Population pharmacokinetic modeling of vedolizumab for graft‐versus‐host disease prophylaxis in adults with allogeneic hematopoietic stem cell transplant

Abstract We aimed to characterize the population pharmacokinetics (PK) of vedolizumab for acute graft‐versus‐host disease prophylaxis in adults undergoing allogeneic hematopoietic stem cell transplantation (allo‐HSCT) and assess potential clinically relevant covariates. Dosing, patient characteristics, and PK from a phase 1b, open‐label, dose‐finding study of vedolizumab 75 mg initial dose escalated to 300 mg and a phase 3 study of vedolizumab 300 mg in patients receiving allo‐HSCT were analyzed using a two‐compartment population PK model with linear elimination. Covariates included age, race, weight, sex, albumin, lymphocyte count, GvHD type, and concomitant medications. Weight, albumin, and lymphocyte count were time‐varying covariates. Model selection was driven by goodness‐of‐fit criteria, precision of parameter estimates, and visual predictive checks. In 193 patients undergoing allo‐HSCT, vedolizumab PK were well described by a two‐compartment, linear PK model. Using reference covariate values, final parameter estimates (95% confidence intervals [CI]) were: clearance, 0.148 (0.136, 0.162) L/day; central volume of distribution, 3.12 (3.03, 3.21) L; intercompartmental clearance, 0.500 (0.408, 0.612) L/day; and peripheral volume of distribution, 3.95 (3.52, 4.44) L. Weight and albumin were the most important predictors of vedolizumab PK, with clearance decreasing by ≈20% for low body weight/high albumin and increasing by ≈30% for high body weight/low albumin. There was an inverse relationship between vedolizumab clearance and age, but no detectable effect for lymphocyte count or GvHD type. Post hoc analyses did not detect any relationship between vedolizumab PK and concomitant medications. In summary, the covariates studied did not have a clinically meaningful effect on the PK of vedolizumab.


| INTRODUC TI ON
Allogeneic hematopoietic stem cell transplantation (allo-HSCT) is a potentially curative therapeutic option for patients with hematologic malignancies, but is frequently complicated by the emergence of graft-versus-host disease (GvHD). 1,2[3] In acute GvHD, tissue damage due to the allo-HSCT conditioning regimen leads to the release of damage-associated molecular patterns (DAMPs) and pathogen-associated molecular patterns (PAMPs) that cause activation of host antigen-presenting cells. 4Donor T cells then become primed in the secondary lymphoid tissue of the host and subsequently migrate to target organs such as the GI tract where inflammatory cytokines and effector T cells damage the gut mucosa, causing the symptoms of acute GvHD. 5,6The trafficking of T cells to GI mucosa and gut-associated lymphoid tissue is mediated through adhesive interactions between integrin α4β7 and mucosal addressin cell adhesion molecule-1 (MAdCAM-1). 7Additionally, integrin α4β7 has been shown to be significantly increased on the surface of naive and memory T cells in individuals who subsequently developed intestinal acute GvHD (aGvHD) following allo-HSCT, versus those who developed skin aGvHD or no GvHD. 5,8Studies in murine models of aGvHD suggest that preventing T-cell trafficking to gut-associated lymphoid tissue through disrupting interactions between integrin α4β7 and MAdCAM-1 may reduce or prevent the occurrence of aGvHD, [9][10][11] indicating that integrin α4β7 represents a viable target for the treatment or prevention of aGvHD.
Vedolizumab (ENTYVIO®, Takeda Pharmaceuticals, Inc.) is a humanized immunoglobulin G1 (IgG1) monoclonal antibody that selectively targets integrin α4β7 and is used to treat patients with moderate to severe active ulcerative colitis or Crohn's disease. 12,13dolizumab inhibits the interaction between integrin α4β7 and MAdCAM-1, preventing the migration of gut-homing leukocytes into the GI mucosa, 14,15 and provides gut-selective, anti-inflammatory effects without systemic immunosuppression. 168][19] Although a phase 2a clinical trial in a small number of patients with advanced, severe disease failed to demonstrate positive proof of concept for the efficacy of vedolizumab treatment in patients with steroid-refractory GvHD, 17 interest is growing in the use of vedolizumab for prophylaxis of intestinal aGvHD. 15,20,21A phase 1 trial has demonstrated that intravenous (IV) vedolizumab 300 mg was well tolerated when added to tacrolimus and methotrexate for GvHD prophylaxis, and patients receiving vedolizumab at this dose experienced a relatively low incidence of lower intestinal aGvHD. 15Additionally, a phase 3 clinical trial which evaluated the efficacy and safety of vedolizumab for the prophylaxis of intestinal aGvHD in patients undergoing allo-HSCT (VEDO-3035/GRAPHITE/ NCT03657160), found that vedolizumab achieved a statistically significant and clinically meaningful improvement in lower gastrointestinal aGvHD-free survival by Day 180 after allo-HSCT. 17eviously, vedolizumab pharmacokinetics (PK) has been characterized in a population modeling study conducted in adults with ulcerative colitis or Crohn's disease. 22Vedolizumab PK was best described by a two-compartment model with parallel linear and nonlinear clearance pathways. 22Extreme albumin and body weight values were found to have a potential to be clinically important predictors of clearance in patients with inflammatory bowel disease (IBD; effect sizes >25%). 22Patients with IBD and those requiring allo-HSCT for malignant and nonmalignant hematologic conditions are different patient populations and the pathophysiology of IBD and hematologic conditions are heterogeneous.
In this study, a previously developed population PK model 22 was refined to describe the disposition of vedolizumab administered for aGvHD prophylaxis in two clinical studies.We aimed to characterize the PK of IV vedolizumab for aGvHD prophylaxis in adults undergoing allo-HSCT, and to assess the impact of potential clinically relevant covariates on vedolizumab PK in these populations.

| Study design and analysis populations
Data were derived from two previously reported clinical studies: VEDO-1015, a phase 1b, open-label, dose-finding study of vedolizumab (NCT02728895) and VEDO-3035 (GRAPHITE/ NCT03657160); a phase 3 study to evaluate the efficacy of vedolizumab when added to background standard of care aGvHD prophylaxis in patients who received allo-HSCT. 15Briefly, the studies were conducted in patients undergoing allo-HSCT for a haematologic malignancy or myeloproliferative disorder (full inclusion and exclusion criteria have been published previously). 17

| Ethics approval statement
The studies were approved by each site's institutional review board, and all participants provided informed consent in accordance with the principles of the Declaration of Helsinki and the International Council for Harmonization Guidelines for Good Clinical Practice.
Nomenclature used for drug and molecular target conforms to the IUPHAR/BPS Guide to Pharmacology nomenclature classification. 23

| Analytical assay
Serum concentrations of vedolizumab were determined using a sandwich enzyme-linked immunosorbent assay method with a validated range of 0.20-8.0μg/mL, as noted previously. 22The presence of AVA in serum was detected using validated electrochemiluminescence (ECL) assays. 24

| Data management
In this analysis, data from the allo-HSCT patients from the phase 1b study (VEDO-1015) and the phase 3 study (VEDO-3035) who were administered at least one dose of vedolizumab and for whom at least one quantifiable post-dose plasma PK sample was available were included in the PK analysis.Dosing, PK, and covariate data from patients who received vedolizumab 75 and 300 mg 15 were merged and assembled into NONMEM®-ready dataset using R (version 4.1). 25ngitudinal vedolizumab serum concentration measurements were used for the PK analysis.Covariates included in the dataset were age, body weight, sex, serum albumin level, absolute lymphocyte count, evidence of aGvHD (liver, skin, or intestine), and anti-vedolizumab status.Weight, albumin, absolute lymphocyte count, and evidence of aGvHD (at any time during the study and of any grade) were all included as time-varying covariates, with interpolation of values between observed time points performed using a next observation carried backward approach.The remaining covariates were patientlevel or baseline factors.

| Population PK model development
The population PK analysis was performed using NONMEM® software (version 7.5, ICON Plc, Ireland).The first-order conditional estimation with interaction method was used for model development. 26rst, a base model was developed without considering covariate effects.The previously developed population PK model in patients with IBD was used as the starting point for the development of the base model in this analysis, after removal of all covariates. 22e structural form of the previous model was a two-compartment model with zero-order input and parallel linear and nonlinear clearance pathways, parameterized in terms of clearance (CL), central volume of distribution (V c ), intercompartmental clearance (Q), peripheral volume of distribution (V p ), maximum elimination rate (V max ), and concentration at half-maximum elimination rate (K m ), with appropriate random effect distributions. 22Additional structural models were tested, including a standard two-compartment linear PK model.Appropriate models for interindividual variability and residual error were also investigated during base model development.Interindividual variability was described using an exponential random effects model with the assumption of a log-normal PK parameter distribution.For the residual error model, both proportional and combined additive and proportional models were tested.
Inter-occasion variability (IOV) was modeled on CL using an exponential random effect.A new occasion was defined whenever a vedolizumab dose was administered with subsequent vedolizumab PK observations taking place prior to the next administered dose.

The assessment of model suitability and decisions about increasing
or decreasing model complexity were guided both by goodness-offit criteria.Goodness-of-fit criteria included diagnostic scatter plots, successful convergence of the minimization routine, plausibility and precision of parameter estimates, correlation between model parameter estimation errors of <0.95 and Akaike information criterion (AIC).Measures of estimation precision included percentage relative standard errors (%RSEs) and 95% confidence intervals (CI) derived from NONMEM®-reported SEs.
Second, a full covariate model emphasizing parameter estimation rather than stepwise hypothesis testing was implemented, owing to typical problems encountered when using stepwise regression techniques. 27The full covariate model was developed by incorporating the effect of all prespecified covariate-parameter relation- The resulting distributions of relative CL for a given covariate value were visualized using a forest plot.
As a hypothesis-generating exercise, an exploratory post hoc analysis was performed to investigate the potential impact of aGvHD on CL.Exploratory graphics included individual random effects for CL from the final model versus type of aGvHD (skin, intestinal, or liver).Additional post hoc analyses explored the influence of select concomitant medications on vedolizumab PK.
Patients had a mean age of 50.8 years at baseline (Table 1).The final population PK dataset was comprised of data from 193 participants from two studies contributing a total of 2380 vedolizumab observations.Of these, 15 samples (0.6%) were below the limit of quantification (0.2 μg/mL) and were subsequently excluded from the analysis.
In addition, four observations were excluded from the analysis because they were either below the limit of quantification immediately after the end of infusion or were unusually high concentrations.Six patients (3.1%) had evidence of liver GvHD at any time, while 62 (32.1%) had skin GvHD and 17 (8.8%)had intestinal GvHD.Two patients had antibodies against vedolizumab during the vedolizumab dosing period; positive tests for both patients occurred on the first day of dosing.At baseline, mean absolute lymphocyte count was 0.426 K/μL.Mean albumin for patients at baseline was within the normal range, at 3.97 g/dL.

| Base PK model
Initially, a two-compartment model with linear and nonlinear clearance was evaluated on the basis of a previous population PK analysis of vedolizumab in patients with ulcerative colitis and Crohn's disease. 22However, the absence of adequate data at low concen-  S1.During model development, various random effects and residual error structures were tested.The final model included IIV, using exponential variance models, on CL, V c , and V p in a full block matrix: IIV percent coefficient of variation (CV%) was 44.0%, 22.5%, and 39.0%, respectively.The IOV was estimated on CL and the CV% was 21.5%.These estimates of random effects parameters show that clearance was moderately variable, both between individuals and between dosing occasions.The residual error model included only a proportional term (CV% = 15.9).All parameters were well estimated with adequate precision as judged by the CI ranges.
Shrinkage was low for IIV of CL and V c (8.31% and 6.40%, respectively), moderately high for IIV of V p (26.4%), and high for IOV of CL (56.1%).The base model provided an adequate description of the data, aside from one particularly high observed concentration (above 500 μg/mL), the population and individual predicted vedolizumab concentrations matched the observations well.No systematic bias was noted in the conditional-weighted residual or NPDE versus time, population predicted concentration, or time after dose plots.aside from a tendency toward low NPDE values for later times after dose, where data were sparse.Similarly, there was no systematic bias in NPDE versus covariate plots, except for some deviations for low weight and albumin values.
The random effects were mostly symmetric around zero, with positive correlations observed between random effects for CL and V c and for V c and V p , and little to no correlation between random effects for CL and V p (Figure S1).There were some obvious trends evident in the relationships between interindividual random effects and continuous covariates: between body weight and CL, V c , and V p , and between albumin and CL.There were also some relationships observed between interindividual random effects and categorical covariates (study, sex, and race), although these are potentially explained by differences in body weight in these groups.There were no apparent trends between inter-occasion random effects on CL and the time-varying covariates, aside from a minor deviation for TA B L E 1 Covariates of study participants contributing vedolizumab observations to the pharmacokinetic dataset.for CL was estimated at 17.9 CV% (Table 3).Most structural model parameters, covariate effects and random effect terms were estimated with reasonable precision based on %RSE and 95% CI results.

Standard goodness-of-fit plots, including model predictions versus
observations, residuals versus model predictions and residuals versus time, revealed that the final model was generally consistent with the observed data and minimal systematic bias remained (Figure 2 and Figure S7).weight (Figure 3).Albumin was observed to be the most influential covariate and CL was observed to increase with decreasing albumin concentration.CL increased in line with increasing body weight, based on the allometrically scaled CL model (Figure 3).Following post hoc exploratory analysis, no clear relationships were evident between CL and type of aGvHD or concomitant medications after accounting for covariates included in the final population PK model (Figure S8).

| Final PK model evaluation
The VPC demonstrated good alignment between model-simulated data and observed data (Figure 4), with the median, 5th and 95th percentiles of the observed data largely falling within the respective 95% prediction intervals of the simulated data.(Figure S9), nor over the ranges of the covariates included in the final model (Figure S9).

| DISCUSS ION
This study is the first to assess population PK of vedolizumab in adult patients undergoing allo-HSCT for GvHD prophylaxis and assess the effect of covariates on vedolizumab PK in this patient population.This analysis was based on pooled data from 24 patients from a previously reported phase 1b study, in which patients undergoing allo-HSCT received vedolizumab 75 or 300 mg plus tacrolimus and methotrexate for GvHD prophylaxis 15 and 169 patients from a phase 3 study to evaluate the efficacy of IV vedolizumab 300 mg when added to background GvHD prophylaxis in patients who received allo-HSCT.
Using these data, vedolizumab PK were well described using a two-compartment model with first-order elimination, and the characteristics of vedolizumab PK were generally consistent with other therapeutic monoclonal antibodies of the IgG1 subclass. 28ditionally, PK disposition and identified significant covariates were consistent with the previous population PK analysis of vedolizumab performed in patients with IBD. 22In the evaluation of patients with IBD, the structural PK vedolizumab model was two-compartment with zero-order input and parallel linear and nonlinear clearance pathways, and was largely consistent with observations from this study, though only linear clearance was selected for use in the current final model. 22In patients with IBD, vedolizumab exposures decreased in a biexponential fashion after reaching peak concentrations, with a biexponential decrease observed down to approximately 1-10 ng/mL, after which concentrations appeared to decrease in a nonlinear fashion.The contribution of nonlinear elimination in IBD patients was, however, minimal to the total clearance of vedolizumab at the approved therapeutic dosage of 300 mg. 22In this analysis, vedolizumab exposure was also observed to decrease biexponentially, though at- The clinical implications of the findings of this study for vedolizumab in a population undergoing allo-HSCT must be considered.A shaded region from −25% to +25% of the reference value has been added to provide visual context for each covariate effect.This shaded region should not be used to make statistical inferences about the magnitude of the covariate effects.

Median [95% CI]
F I G U R E 4 Final model: Prediction-corrected visual predictive check of vedolizumab concentration versus time after most recent dose.Black lines represent the median (solid), 5th and 95th percentiles (dashed) of the observed data.Blue and grayshaded regions represent the 95% confidence intervals of the corresponding (i.e., 5th, 50th, and 95th) percentiles.
ships into the base model, while avoiding simultaneous inclusion of moderate to highly collinear/correlated covariates (absolute values of correlation coefficients >.30).Determination of the prespecified covariates was based on clinical interest and previous knowledge of the clinical pharmacological profile of vedolizumab.22These covariate-parameter relationships included the effect of albumin on CL, the effect of body weight on CL, V c , Q, and V p , and the effect of absolute lymphocyte count on CL.The full covariate model was constructed by simultaneously including the effects of these prespecified covariates in the base model (EquationS1).The impact of body weight on vedolizumab PK was described using an allometric model with a reference weight of 75 kg (estimated and fixed allometric exponents were investigated).Other continuous covariates were included in the model as power functions normalized by a reference value.Categorical covariates were included in the model as power functions, with a dichotomous (0, 1) power coefficient serving as an on-off switch for each effect.Inferences regarding the clinical relevance of each covariate from the full model were based on both the resulting parameter estimate and estimate precision.Uncertainty distributions of estimates for CL and the covariate parameters were constructed by simulating from a multivariate normal distribution derived from the NONMEM® covariance matrix of estimates.The median simulated CL estimate for the reference individual was designated as the "reference CL." Distributions of the reference CL adjusted for a given covariate value were constructed on the basis of the median and 90% CI from the uncertainty distribution of the covariate-parameter estimates.The covariates were evaluated in a univariate fashion (i.e., covariates were fixed at their reference value except when perturbed).All derived CL values were normalized to the reference CL to reflect the relative change in CL versus the reference individual.
Standard goodness-of-fit diagnostic plots, simulation-based visual predictive checks (VPC), and normalized prediction distribution errors (NPDE) were used to assess the predictive performance of the model.For the VPC, 500 Monte Carlo simulation replicates of the original dataset were generated using the final population PK model.Observed vedolizumab concentrations were plotted versus time after the first dose, including the observed median and 5th and 95th percentiles.Observed plots were then overlaid with the 95% CI of the median, 5th and 95th percentiles of the simulated data for comparison.NPDE were plotted against population-predicted vedolizumab concentrations and time after the first dose to check for evidence of model mis-specification and were also plotted against covariates included in the model to gauge the presence of any bias.Histograms of the NPDE were then plotted to evaluate their distribution, which should be normal with a variance of 1 and a mean of 0.
trations (only three patients received the 75 mg dose) precluded adequate estimation of the nonlinear clearance parameters (V max and K m and model development continued with only linear clearance).A two-compartment, linear population PK model continued to provide the best and most parsimonious fit to the data.The base model was implemented in NONMEM® with the ADVAN3 TRANS4 subroutine, and was parameterized in terms of CL, V c , Q, and V p .The base model also included fixed allometric weight effects on all clearance and volume parameters.Both models resulted in similar AIC values and goodness-of-fit diagnostics, indicating that the data were adequately described by the more parsimonious two-compartment model with linear CL only (selected base model), compared with the two-compartmental model with linear and nonlinear clearance.Composite plots of vedolizumab concentrations over time after the first dose and time after the most recent dose are shown in Figure 1.Peak vedolizumab concentrations appear to increase with time over the first several occasions before reaching steady state.The parameter estimates for the base PK model are shown in Table low albumin values, where data are sparse.Distributions of the IOV did not appear to vary over time.Plots of individual and population predictions of vedolizumab concentrations versus time matched the observed data well.The prediction-corrected VPC for vedolizumab versus time after last dose showed a reasonable ability to predict the observed data overall.Diagnostic plots for the final base PK model are shown in Figures S2-S5.3.2.2 | Final PK model (full model) The final covariate model expanded the base model to include covariate effects.The final model included the effects of body weight on all clearance and volume terms, with estimated exponents for CL and V c and fixed allometric exponents for Q (0.75) and V p (1).Based on plots of random effects, the effect of albumin on CL was also included in the model.Effects of study, sex, and race were not included in the final model due to their association with bodyweight.The final model also included effects of age, lymphocyte count, and GvHD (liver, skin, and intestine) on CL, although there were no clear relationships observed between these covariates and random effects.The effect of antidrug antibodies on vedolizumab PK was not evaluated due to the limited number of subjects with positive ADA tests.Plots of individual model fits for the final vedolizumab PK model are provided in Figure S6.In the final population PK model, parameter estimates were 0.148 L/day (CL), 3.12 L (V c ), 3.95 L (V p ), and 0.500 L/day (Q), using typical reference covariates for body weight (75 kg), albumin (4 g/ dL), and absolute lymphocyte count (0.1 K/μL; Table2, Equation S2).Final estimates of unexplained variability in CL, V c , and V p were 29.4 percentage coefficient of variation (CV%), 18.1 CV%, and 44.2 CV%, respectively.Shrinkage values for the CL and V c , interindividual random effects were small with estimates of 12.3%, and 8.82%, respectively (estimates based on the standard deviation scale).IOV The relative changes in CL for a reference individual with various covariate values are illustrated in Figure 3. Median and 90% CI values for the reference individual reflect the uncertainty around the population estimate of CL for an individual with reference covariate values.For evaluation of covariates, the point marker and interval lines reflect uncertainty in the univariate covariate effect estimate and the change from the reference CL (i.e., the median posterior CL estimate for the reference individual).Covariates that were identified as having potentially clinically relevant effects on CL included extreme values of albumin and body F I G U R E 1 Composite plots of vedolizumab concentrations versus time after first dose in (A) VEDO-1015 and (B) VEDO-3035 (GRAPHITE).Occasion: A new occasion was defined whenever a vedolizumab dose was administered with subsequent vedolizumab PK observations taking place prior to the next administered dose.
modeling of the nonlinear clearance component caused the structural model to collapse to a linear clearance (first-order)only model.The parameters V max and K m were therefore not able to be estimated well in the absence of data at low concentrations.It is possible that nonlinear clearance of vedolizumab occurs in an allo-HSCT population; however, these data were unable to support this, potentially because the majority of observed vedolizumab concentrations were >10 ng/mL, and limited PK data were available for the lower vedolizumab dose of 75 mg owing to the small number of patients in this group (n = 3).Importantly, the final PK model adequately described the observed data as judged by VPCs, NPDEs, and standard goodness-of-fit plots.Covariates investigated in this study were chosen on the basis of previous modeling in patients with IBD and factors of specific clinical interest in populations undergoing allo-HSCT (such as lymphocytes).22Linear CL was found to be 0.159 L/day in patients with ulcerative colitis and 0.155 L/day for those with Crohn's disease, similar to that observed in the current population undergoing allo-HSCT (0.148 L/day).22In addition, relative similarities were also observed among the other PK parameters between the previous and current analyses, such as central volume of distribution with values of 3.19 and 3.12 L estimated for V c 22 ; however, minor differences were observed between PK parameters of V p and Q, which were higher in the GvHD model compared to the IBD population with values of 1.65 and 3.95 L estimated for V p ; and 0.12 and 0.5 L/day estimated for Q, respectively.Larger volumes of vedolizumab partitioning to the peripheral compartments could potentially account for the higher steady-state concentrations observed in the GvHD population compared to the IBD population.This study identified that extreme values of albumin and body weight may have impact on CL as also observed in the earlier study of patients with IBD.22 Final model parameter estimates and measures of estimation precision (95% confidence interval [CI]) were used to infer the clinical relevance of covariate effects based on a reference covariate group.In line with the previous population PK model in ulcerative colitis and Crohn's disease, a covariate effect size on vedolizumab clearance >25% was considered clinically meaningful.22

F I G U R E 2
Observed versus population-predicted (A) and individual-predicted (B) vedolizumab concentrations.Observed values are indicated by black circles, the dashed blue line is a LOESSs smooth, and the line of identity (solid gray) is included as a reference.Notably, when all other conditions were held constant at the reference values (53-year-old patient weighing 75 kg with no GvHD in the liver, skin, or intestine, an albumin value of 40 g/L, and a lymphocyte count of 0.1 K/μL), vedolizumab AUC ss increased by about 20% for subjects with low body weight (10th percentile: 51 kg) or high albumin (90th percentile: 46 g/L), and AUC ss decreased by about 25% for subjects with high body weight or low albumin.At observed body weights ≥105 kg (90th percentile) or albumin concentrations 36 g/L (25th percentile), these covariates had the potential to be clinically important, based on effect sizes and 90% CIs resulting in a change in CL of >25%.Absolute lymphocyte count, age, and the presence and type of aGvHD had no clinically meaningful impact on CL, though this was anticipated because there were few participants representing each aGvHD subtype, which likely precluded full evaluation of this covariate effect.Additionally, only two subjects in the phase 3 study had antibodies against vedolizumab during the vedolizumab dosing period, though this may primarily be due to the administration of myeloablative therapy and immunosuppressant treatment following the first dose of vedolizumab.The scarcity of patients with AVA during the trial meant that the presence of AVA could not be evaluated as a covariate in this study.The findings of this study suggest that albumin level and body weight are statistically significant covariates that affect clearance, as observed in patients with IBD, though these parameters are likely to be relevant only at extreme values.22The post hoc exploration of effects of concomitant medications on vedolizumab PK revealed no relationships between vedolizumab disposition and methotrexate, tacrolimus, cyclosporine, and ursodeoxycholic acid.Importantly, this study has illustrated F I G U R E 3 Covariate forest plot for relative vedolizumab clearance (CL) based on the final model.The colored circles represent the median and the solid horizontal lines represent the 95% confidence interval.Relative CL is shown for the 5th, 10th, 25th, 50th, 75th, 90th, and 95th percentiles of the observed values of body weight, age, albumin and lymphocyte count in the analysis dataset.

Covariates Number of patients, n (%) Study VEDO-1015 (n = 24) Study VEDO-3035 (n = 169) All data (n = 193)
Abbreviations: GvHD, graft-versus-host disease; SD, standard deviation.a Represents a summary of counts based on the type of GvHD reported at any time during the study.An occurrence of GvHD included any grade of GvHD (Grades I-IV), and a patient could exhibit one or more types of GvHD.